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INTRODUCTION 


■"This.! report ts a summary of work coinple ted. under NASA grant NCR 11- 
002-179 entitled "Determination of the Effects of Nozzle Nonlinearities 
Upon the Nonlinear Stability of Liquid Propellant Rocket Motors". Research 
activities supported by this grant were begun in -August 1973^ and satis- 
factory progress has been made toward meeting the research objectives dur- 
ing the first year of effort. Before giving a description of this progress, 
the motivations and objectives of this research project will be briefly 
reviewed. 

’ :\ ^Various aerospace propulsion devices, such as liquid and solid pro- - , 

pellant rocket motors and air breathing Jet engines, are often subject 
to combustion instabilities 'which are detrimental to the performance and 
safety of operation of these devices. .In order to design stable engines.,. ' 
capabilities for a. priori determination of the linear and -nonlinear 
characteristics of the instability and the range of operating conditions 
for which these engines are dynamically stable must be acquired. In 
order to perform such an analysis, the behavior of the exhaust nozzle 
under oscillatory flow conditions must be understood* In particular, it 
is necessary to know how a wave generated in the combustion chamber is 
partially transmitted and partially reflected at the nozzle entrance. 

This information is usually expressed as a boundary condition (usually 
referred to as a Nozzle Admittance Relation) that must be satisfied at the 
nozzle entrance. 

Before such a boundary condition can be derived, the nature of the 

wave motion' inside the nozzle must be investigated. The behavior of 

oscillations in a converging-diverging supercritical nozzle was first 

treated by Tsien" 1 * who considered the case in which the oscillation of 

p q 

the incoming flow is one-dimensional and isothermal. Crocco extended 
Tsien's work to cover the more general cases of non-isothermal one- and 
three-dimensional oscillations. The analyses of Tsien and Crocco are 
both restricted to small- amplitude (i.e., linear) oscillations. More 
recently, a nonlinear nozzle theory has been developed by Zinn and Crocco^ 
who. extended the previous linear theories to the investigation of the 



behavior of finite -amplitude waves . 

In recent studies (supported under NASA grant NGL 11-002-083) con- 
ducted by Zinn, Powell, and Lores, theories were developed which describe 

17- g Q 

the nonlinear behavior of longitudinal' 5 and transverse^ 5 instabili- 
ties in liquid -propellant rocket chambers with quasi-steady nozzles. These 
theories have now been extended to situations in which the instabilities 

are three-dimensional and the rocket combustors are attached to conven- 
11 

tional nozzles . All of these theories have successfully predicted the 
transient behavior, nonlinear waveforms, and limit-cycle amplitudes of 
longitudinal and tangential instabilities in unstable motors . 

A new nonlinear nozzle theory is needed for the following reasons. 

5 6 

First, the nonlinear analysis of Zinn ? is mathematically complicated 
and requires considerable computer time. For this reason, Zinn’s analysis 
has never been used to perform actual computations of the wave structure 
in the nozzle or the nonlinear nozzle response. Secondly, the nonlinear 
nozzle admittance relation developed by Zirni is not compatible with the 
recently developed nonlinear combustion theories (see References 7 through 
11) . Consequently, a linear nozzle boundary condition or short nozzle 
(quasi-steady) assumption had to be used in all of the combustion in- 
stability theories developed to date. With the exception of a few special 
cases, where the amplitude of the instability is assumed to be moderate 
and the mean flow Mach number is small (e.g., see Reference 9) > the use of 
a linear nozzle admittance relation in a nonlinear stability analysis is 
obviously inconsistent. Furthermore, in the case of transverse insta- 
bilities the "linear" nozzle has been known to exez’t a destabilizing 
effect; in these cases it is especially important to know how nonlinearities 
affect the nozzle behavior. 

Tie objective of this research program is to develop a three- 
dimensional, nonlinear nozzle admittance relation to be used as a 
boundary condition in the recently- developed nonlinear combustion insta- 
bility theories .' This objective will be accomplished by performing the 
following four tasks : 

Task 1: Development of the theory 

Task II: Calculation of the nozzle response 

Task III: Application of the nozzle theory to combustion 

instability problems 



Task IV: Preparation of the final technical report 

.During the first six months of this projects considerable progress 
was made toward completing the first -of the above tasks. However , un- 
- foreseen difficulties in. the mathematical formulation of the problem arose 
in December, and most of the first year was needed to complete Task X. 
Once the theory and computer programs were developed, Task IX was . com- 
pleted during the remaining time. A one-year extension of support has 
been granted by NASA to complete Tasks III and IV. A summary of the work 
completed on Tasks ! and XI is given in the remainder of this report.. 


TASK I: DEYELOPTOIT OF THEORY .. 

Derivation of the Nozzle Wave Equation 

5 6 

As in . the Zinn-Crocco.. analysis , finite-amplitude, periodic, oscil- 
lations inside the slowly convergent, subsonic portion of . an ■ axi sym- 
metric nozzle operating in the supercritical range were investigated. . 
The flow in the nozzle was assumed to be adiabatic and inviscid and to 
have no body forces or chemical reactions. The fluid was also assumed to 
be calorieally perfect. ' : . • .. .. .. 

The nondimens ional equations describing the gas motion in the nozzle 
were written in the following form: 


If + V • (pV) = 0 


( 1 ) 


3V 

at 


+ |v (V • V). + (Vx V)xV + ^Vp =0 


( 2 ) 


+ v * VS = 0 

at 


( 3 ) 


S = — top - top + constant 


(V) 



where y is the specific heat ratio; V, p, p, and S are the dimension- 
.less velocity, pressure, density and entropy respectively and t is the 
dimensionless time. 

It was also assumed that the nozzle flow is isentropic and irro- 
tational.. Under these conditions the energy equation (i.e.. Equation 
( 3 )) is no longer needed, the state equation (i.e.. Equation (4)) reduces 
to the isentropic flow relation, p = p , and a velocity potential exists 
such that « V. Hie continuity and momentum equations were combined, 
with the aid of the isentropic relation, to yield the following equation 
which describes the behavior of the velocity potential: 

P 2 

V $ - + (y ** 1) V § ( 5 ) 

+ ,, !Y„ ,. r . .J r . (v$*v$) v 2 § + ~ vf -v(v§ -v$) 

Hois equation is consistent with the wave equation used in the second-order 
nonlinear combustion instability theory developed by Powell, Zinn, and 
Lores (see References 7 and 10) . 

In the nonlinear combustion instability theories developed by Powell, 
and Zinn, each variable was expressed as the sum of a space -dependent 
steady state quantity and a time- and space -dependent perturbation quantity 
In order to obtain a nozzle admittance relation compatible with these theo- 
ries, the velocity potential was expressed as follows: 

§ = i + ? ' (6) 

where the prime denotes the perturbation quantity and the bar denotes 
the steady- state quantity. Using the relation V® = V, Equation (6) was 
substituted into Equation ( 5 ) to obtain the following wave equation for the 
nozzle: 

[1 - V 2 ] v 2 ®' - = 2 v- [v®^ + | V(V-V$')] 


+ (7) 
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(y - 1) (v-v) + V* ']. + | vCv 2 ) * 


7? 


2 V$' • [v^ + | v(v*v#')] + | V-vCvf'-Vf') 


+ (y - l)v 2 $' + v-v§'j + Y ~ — (v*v) (vs'-vO 


1 r. A / a . ?/■' 


+ -| v _ “ JL v 2 § '(v§ 7 -v $ ') + | v§'*v(v?'*v? / )} 


.Before proceeding, 'with the analysis, a coordinate system, appro- 
priate, for the introduction of the boundary condition at the nozzle 
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vails', was chosen. * 'Following the approach used by Zinn and Crocco 5 
for an axi- symmetric nozzle, the axial variable z was replaced by the 
steady-state potential function cp, and the radial variable r was re- 
placed by the steady-state stream function . Hie potential and stream 
functions are defined by: 


d\jj 

rpu = ^ 



( 8 ) 


where 6s and 6n respectively represent elementary (non-dimensional) 
lengths in the directions of the unperturbed streamlines and of their 
normals on the meridional planes (see Figure l) and u is the steady- 
state velocity. A third independent variable, 8, measures the azimuthal 
variation. In the new coordinate system, the perturbation velocity is 
expressed in terms of its components along the coordinate directions as: 


„/ / , t . r 

V - u e + v e . + w e ^ 
— • cp ij; -* 8 


( 9 ) 


where the e f s are unit vectors. 

The transformation of Equation (7) to (qr,4,8) coordinates was greatly 
simplified by assuming that the steady-state flow is one -dimensional, which 
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Figure 1. Coordinate System used for the Solution of the Oscillatory Nozzle Flow. 





is a good approximation- for slowly convergent nozzles. Under these con- 
ditions the dependence o£ p and u on f and 9 can be neglected, so that 
they are considered t.o be practically uniform on each surface <p ~ constant. 
Also the angle of obliquity of the stream- lines to the axis of symmetry 
is sufficiently small so that its cosine is practically 1 and the. element 
of normal fin along the surface cp - constant can be identified with dr. 

Hence the first of Equations (8) was integrated to obtain: 

r 2 = “ • - * (10) 

P u 

In addition the mean ..flow . velocity (vector appearing in Equation (7) is . 
given by: 


V = u (cp) e 


cp 


( 11 ) 


With the aid of Equations (10) and (ll) and the expressions for the 
Laplacian, divergence,, and gradient in a (<p,ijr,9) coordinate system, 
Equation ( 7 ) was transformed to the following equation: 

f l (!p) ■ f 2 (tp) ^ + f 3 (q,) [ 2 U*n + **) + i *9e] (12) 


- 2 Kt + K - h ?, tt 

u 


■ 2 + f hi + j ? e" ? 9't 

+ ( Y + X)u 2 ^ + 2 pu **' 


+ f 1 + f (g) ( 5 ') ^ 

2 ^ 0 Gcp 


+ f 6 (<p) t (i') 2 + f 6 (,) jT ($0 ') 2 + ( Y - 1 ) 
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- V«p). »«& + < Y - d| [2 (n' n * *;) 


t 

s e 9 ] K + (y " «" 4 2 (+*♦* + *♦) 


2+ e J "cp 




f 2 (cp) 


-2 -2 
c - u 


„ -2 
1 du 

-2 dcp 


f 3 ( <p) 


f 4 (<p) = 


-2 dcp 
2c 


;<*> = I t 


•, “2~i 2 

- 1 u du 

2 -2J. dcp 

c 


f 6 (q>) = -E- [1 - (2 - y) 


- 2 n -2 
u 1 du 


In Equations (13) c is the steady-state sonic velocity given by: 

-2 v - 1 -2 

c » 1 - — u 


In deriving Equation (12) the third-order terms in Equation ( 7 ) (i.e., 
the last two terms on the right-hand side) have been neglected, thus 
Equation (12) is correct to second order. 

Application of the Cfalerkin Method 

... >, , , ■ 

The equations obtained by the above procedure have no known closed- 
form mathematical solutions. Consequently, it is necessary to resort to 



the use of either numerical solution techniques or approximate analytical 
techniques. Since the numerical solution techniques generally require 
excessive computer time, the latter approach was used. In the nonlinear 
combustion instability theories developed by Powell and ZInn (see Refer- 
ences 7 Jl) the governing equations were solved by means of an approxi- 
mate solution technique known as the Galerkin Method* which is a special 
case of the Method of Weighted Residuals > . In these investigations 

it was shown that the Galerkin Method could be successfully applied in the 
solution of nonlinear .combustion instability problems; its application 
was straightforward and it required relatively little computation time. 

Thus the Galerkin Method was also used in the present analysis to de- 
termine the nonlinear nozzle admittance relation. 

In. order to employ the Galerkin Method in the solution of the wave 
equation (i.e., Equation’ (12) * it was first necessary to express the ve- 
locity potential, $ f , as an approximating series expansion. The struc- 
ture of this series expansion was guided by. the experience gained in 
the nonlinear nozzle admittance studies performed by Zinn and Crocco 
(see Reference 5) as well as in the nonlinear combustion instability 
analyses of Powell and Zinn (see Reference 10) . Thus the velocity po- 
tential was expressed as follows: 

M N 

= I X cos m9 

m=0 n=l 

In Equation (15) * the functions A (cp) are unknown complex functions 

of the axial variable cp. The 9- and ^-dependent eigenfunctions were de- 

termined from the first-order (i.e-, linear) solutions by Zinn . In these 

functions m is the transverse mode number, n is the radial mode number, J 

m 

is a Bessel function of order m, is the value of. the steady-state stream 

w 

function evaluated at the nozzle wall, and S is a root of the equation 

mn . 

d-J (x) /dx = 0 . The expansions given above describe standing wave motion; 

they can be easily modified to describe spinning wave motion. In the time- 

dependence, (ju is the fundamental frequency which must be specified and the 

integer k gives the frequency of the higher harmonics . The values of k 

mn mn 

for the various modes appearing in Equation (15) must be determined from 


h (- 

L mn \\ 


qpif] 

m L mn \f j 


'w 


ik u)t\ 
rnn J 


(15) 



the results of the nonlinear combustion instability analysis of Powell 
and Zinn^" . For example it was found that, due to nonlinear coupling be- 
tween modes, the second tangential (m = 2, n = l) and first radial (m - 0, 
n 1) modes oscillated with twice the frequency of the first tangential 
(in - 1, n = l) mode.- Thus in Equation (15) k -j 1 ~ 1 for the first tan- 
gential mode and k^ = 2 for the second tangential and the first radial 
modes . The amplitudes and phases of the various modes depend on the 
axial location (i.e., cp) in the nozzle through the unknown functions 

A (cp) . 
run 

In order to simplify the algebra involved in the application of the 
Galerkin Method, the approximating series expansion for § is written as 
a single summation as follows : 

r-, ik u)t 

i' = l A p (<p) 0p(0) Yp«i) e p (16) 

P=1 


where to each value of the index p, there corresponds the mode numbers 

m(p) and n(p) , which determine the value of k . In Eq. (16) © (9) and 

3? P 

¥ (^f) are the 8 -and ^-dependent functions while 35 is the number of terms 
In the series expansion. In the present analysis, a three-term expan- 
sion consisting of the first tangential ( p = 1; m = 1, n = 1) , second tan- 
gential (p = 2; m = 2 ; n = 1) and first radial (p = 35 m = 0 , n = 1) modes 
was used, but the theory is applicable to any number of modes. 

In order to obtain the solution, the unknown cp-dependent functions, 

(cp) , were determined by the Galerkin Method as follows. The assumed 
series expansion for the velocity potential (i.e., Eq. (16)) was sub- 
stituted Into the wave equation to form the residual, E (TO . In the 
event that this residual is identically zero, the assumed solution is 
an exact solution. The residual, therefore, represents the error in- 
curred by using the approximate solutions given by Eq. (16) . The Galer- 
kin Method determines the amplitudes A (cp) that minimizes the residual 

e(T') . P 

Applying the. Galerkin Method, the residual E(T0 was required to 
satisfy the following Galerkin orthogonality conditions': 
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The weighting functions T.(t), ®.(6) and Y . ( t) correspond to the terms 
that appear in the assumed series expansion . The temporal weighting 
^function, T (t) , is thev complex conjugate -Qf.. the- assumed time dependence-,. 


;nur 


-ik cut 

T.(t) = e ^ 


( 18 ) 


The azimuthal weighting functions, ®.(9) , .are given by 

0 


©.'(0) = cos m9 
J 


(19) 


while the radial -weighting functions, Y.(t) ? are given by 

0 


V* } = 


(£ 0 ) 


The time integration ’.is performed over one period of oscillation, 

Ott s 

T = while , the spatial integration is pei'formed over any surface 
of cp = constant in the nozzle (in Eq.' (l7) dS indicates an incremental 
area on this surface) . 

Evaluating the spatial and temporal integrals in Eq. (l?) yields 
the following system of N nonlinear, second order, coupled, complex 
ordinary differential equations to be solved for the complex amplitude 
functions, A (cp) : 


M 


I fe 

P-1 • 


d 2 A (cp) 

p - 

dco 


d/\ (cp) 
+ CL — £ 


dcp 


+ C 3 A p (cp) 


N 

V 

+ L 


N 


I fe 


d 2 A (cp) d 2 A (cp) dA (cp) 

— V V*> + d 2 — lr- 


dcp 


p=l q=l 

<3A (T) dA.(to) 


dco 


dA„(cp) 


+ D, 


H — — + D, A (cp) _ 

3 dcp . dcp 4 p v dcp 


+ D_A 

5 p 


(cp) A (cp) } 


+ Q — 0 , 


j - 1,2, . * . N • * 


(21) 



In the above equations, Q represents the additional nonlinear terms 
that arise when a complex solution (i.e. Eq. (1 6 )) is used to solve the 
nonlinear wove equation (i.e. Eq. (12)) . These terms are similar in 
form to the nonlinear terms shown, but they. involve the complex con- 
.j ugates of the . amplitude functions . The procedure for deriving these 
terms is given in Appendix B of Ref. 11. The coefficients C R and are 
functions of the axial variable ip as well as the indices j ,p and q. 
Analytical expressions for these coefficients contain integrals involv- 
ing trigonometric and Bessel functions . In the absence of closed- 
form expressions for the integrals of Bessel functions, these integrals 
were computed numerically . 

•As a check on the above analysis, a single mode series consisting 
of the first tangential mode was used in deriving Eq. (21) . For this 
case, all the coefficients - of the nonlinear terms vanish and the re- 
sulting linear equation is: 


- 2,-2 - 2 s 

u (c - u } 


d"A -2 f 1 

, 2 “ U Us 

dcp c 


du 2 

dcp 


+ 2imj 


dA 

dcp 


( 22 ) 


r n —2 y 

1 " 2T pUG " ~ 


-2 -2 

- 1 . u du ^ 2 

— io) i — + & 

2 -2 dcp 


'} A (<f) 


= 0 


Q 

which is identical to Crocco and Sirignano's equation for the isentropic 
and irrotational case. 


Dominance of the IT Made 

The well known fact that most transverse instabilities behave like 

the first tangential (It) mode was used to further simplify Eq. (21) . 

H 

Based on the results of the recent combustion instability theory, it 
was assumed that the amplitude of the IT mode was considerably larger 
than the amplitudes of the remaining modes in the series solution. 
Through an order of magnitude analysis, correct to the second order, 

Eq. (21) reduced to the following system of equations: 
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2 

A -2 

-2,-2 -2x 1 -2 r 1 du , 0 . | 

u (c - u ) — g - - u Ll 2 + 2ira J 

dcp c 


, dA.. 

! 1 

dcp 


[- 


3 1 2 y -'1 


-2 -2 
u du 


et p uc - -*- 5 - -=•—- + ™ 


W 


5 2 d< P 


V«p) 


= 0 


(23a) 


- 2,-2 

u (c 


„ d 2 A 
- 2 ' _£ 


- u ) 


dcp" 


-2 r 1 du 2 . ' 

u ”3 — + 2ik u) 

L-2 dcp p . 


dA 

dcp 


{■ - 


p — 2 
pUC 


w 


v - 1 U 2 du 2 '. , V 2,21 A ( ^ 

xk UJ “7T -r — + k 0) A (coj 

2 p ^2 dcp p J P ‘ 


2 

a A - 


d 2 A 1 dA 1 


“ ?1 W»v) 2 A 1 “ D 2^ ?p ) 2 dcp 

dcp dcp 


/dA-.s dA^ 2 

~ D 3 (9»p) \“dcp“j " d 4 ^>p) “Hep" A 1 ~ a i 


- Q = 0, , 

"T> 5 


p = 2 , 3 , vi<N- 


The above equations can be written concisely as follows: 


(23b) 


H p (cp) 


d 2 A (cp) 

E_ 

, 2 
dcp 


dA (cp) 

+ M (cp) — " P — — + N (cp) A (cp) = I (co) 
P ^ dcp v ^ P P ' 


( 24 ) 


P = !? 2 ? W , 

where (cp) = 0 . 

Xt can be seen that the above equations are decoupled with respect 
to the XT mode; that is, the solution for A^ can be obtained indepen- 
dently of the amplitudes of the other modes. Thus, to second order, the 
nonlinearities of the problem do not affect the IT mode. On the other 
hard, the nonlinearities influence the amplitudes of the higher modes 
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(i . e - > A 0 , A.,...) "by means of the inhomogeneous terms in the equations 
for the other modes. 


Homogeneous and Particular Solutions • 

• Equation (24) is. a second order, linear ordinary differential 
equation and its general solution is a combination of the homogeneous 
solution that satisfies the homogeneous part of Eq. (24), i.e.. 



d 2 A (h) 

. . , £ , . 
. 2 
dcp 


dA 


(H) 


+ M 


P dcp 



= 0 


(25) 


and the particular solution that satisfies Eq. (24) . The general so- 
lution can be written in the following form: 


A (cp) = q + K^ h) t 


where A and A ^ are two independent solutions of Eq. (25) ? 3C, and 
P P (i) X 

K are arbitrary constants, and A v 1 is a particular solution of the in- 
2 P 

homogeneous equation. 

Examination of the coefficients of Eq. (24) show that this equation 
has the following singular points : 


u = 0 

U SS C 


= (: 


2 ^|2 

,y 4- 1 / 


- C 


throat 


U = cq 


For a supercritical nozzle with a finite area entrance, only the singu- 
larity at the throat is of concern to us. Assuming that the singularity 

of the solution appears in , the condition requiring the regularity 

P 

of the solution at the throat can be expressed by requiring = 0. 
Consequently, the required solution of Eq. (24) is of the form 


A (to) = KjA^ 1 (cp) + A ^4 (tp) 


(26) 
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Derivation, of Admittance Relations 

..... - M 1— -■■■ — , — ■ ... M , ■ 

Using the above result, a nonlinear admittance relation to be used . 
as a boundary condition in nonlinear combustion instability analyses can 
be derived.- Denoting, the terms of Eq. (16) by 

ik ojt 

- A (cd)© (0)1 (f; e P . , ' (2?) 

P P . P p 


talcing partial derivatives with respect to z and t, and using Eq. (26) 
gives 

(i) 


a§' 


“3T? “ u ®^( e ) ^,(t) e P 


ik iDt dA 




dcp 


ik uat dA 

= K-. u 0 (e) 1 (f) e p . 

1 P p . dcp 


(h) 


(28) 


a?' 


ik cut 


- ik (ii 0 ( 0 ) 1 ( )j?) e P A 


(i) 


5t 


P P 


ik u)t / . v 

= K n ik cd 0 ( 0) 1 (\|f) e p A 1 ' 

1 p p w p . p 


. (29) 


Eliminating between Eqs.(28) and (29) and defining 


C - 


af h) /dcp 


(30) 


r = 


r-i ^ h 

— _ s A -1-— 

l p dcp 


- A 


/. v dA^ 
(A) P 


dco 


■] 


( 31 ) 


iuC 

Y — 


p 


VK CD 

P 


(32) 


yi elds 


M' 

+ y y 

A rr. ri nr. 


p ik art 

- uc © (6) 1 (ft) e P T 
P • P Y * ■ P 


( 33 ) 


p = 1,2, — N 


16 


Equation ( 33 ) is the non! '.near nozzle admittance relations to. he 
used as the boundary condition at the nozzle entrance in nonlinear com- 
bustion instability analyses. The right -hand- side of this equation 
arises from the nonlinear terms in the nozzle wave equation. The quan- 
tities Y and T are. respectively the linear and nonlinear admittance 
P P th 

coefficients for the p mode. The nonlinear admittance, r , represents 

-fcr 

the effect of nozzle nonlinearities upon the nozzle admittance and it is 
identically zero when nonlinearities are not present. 

It can easily be shown that Eq. (33) can be written in terms of the 
pressure, and axial velocity perturbations as: 


U - Y P = - uc 2 f p - 1*2*... N (3*0 

P P P P , 


where U and P^ are the amplitudes of the axial velocity and pressure 
perturbations respectively as given by: 


P 


> P (cp) ® (e) v p (*) e p 

P=1 


( 35 ) 


u 


/ 


Y U («)) 0 p (e) p 

P=1 


( 36 ) 


Equation (34) is equivalent to Eq. ( 33 ) to second order only when the 
Mach, number at the nozzle entrance, u , is small. 

In order to use the admittance relation (Eq. ( 33 ) or (34)) in the 
combustion instability theories, the admittance coefficients Y (or C ) 

Jr lr 

and r must be determined for a given nozzle. The equations governing 

these quantities are readily derived from Eq. (24) using the definitions 

for C (i.e., Eqs.(30) and (31)). The resulting equations are: 

P 


K 

H —2. - _ M Q - W 
P dip p p p 



(37) 



IT 


ar 


p &p 


h c 

p p 


+ it 


- i 


p 


25 s 


du 

dcp 
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■ ‘ " p = 1 , 2 , ... N 

i 

TASK II: CALCULATION OF THE NOZZLE RESPONSE 


( 38 ) 


To obtain the nozzle response for any specified nozzle, Eqs. (3?) 
and (38) are solved in the following manner. As pointed out earlier, 
the nonlinear terms vanish for the IT mode (i.e., = 0,1^ « 0) and 

it is only necessary to solve. Eq. (37) to obtain (and hence Xj) at 
the nozzle entrance. Since Eq. (37) does not depend on the higher modes, 
it can be solved independently for Once .has been determined, both 

Eqs. (37) and (38) must be solved for the other modes. In order to do 
tills, the amplitude A^ fcp) must be determined since Eq. (38) depends on 
A^(cp) and its derivatives through I^(cp) . Once (cp) is known, ( 9 ) is 
determined by" numerically integrating Eq. (30) where the constant of in- 
tegration is determined by the specified value of the pressure ampli- 
tude (of the IT mode) at the nozzle entrance. The value of A^ thus 
found is introduced into Eq. (38) which is then solved for r . 

It may be observed that Eq. (37) and (38) have singularities at 
the same points as Eq. (24) . As before, the only singularity of interest 
is the throat. Since Eqs. ( 37 ) and (38) are first order ordinary differ- 
ential equations, the numerical integration of these equations must start 
at some initial point where the initial conditions are known, and termi- 
nate at the nozzle entrance where the admittance coefficients Y and P 

P P 

are needed. Since the equations are singular at the throat, the inte- 
gration is initiated at a point that is located a short distance up- 
stream of the throat. The needed initial conditions are obtained by ex- 
panding the dependent variables in a Taylor series about the throat . 

(cp = 0 ) ; thus , 

C p (cp) = C p (0) + (PC p (0) + ... (39a) 


r (tp) = r (0) + cor'(o) + ... 

p T p p > . 


(39b) 



can be determined by substituting 


The coefficients £ p (°) and Cj(0) 

Eq... (g 9 a) in Eq. (3T)? and taking the limit as cp 0. The results are: 

N (O) 

S (0) " - iHoT 

ir 

-M'(°) C (0) - H'(0) d(0) - »£(0) 

r'( 0 ) = £ £ 2 2 

' P H'(o) 4 - M (0) 


(to) 

(to) 


P = 


... N 


Similarly, F p (0) and l p (o) can be determined by substituting Eq. (39b) 
in Eq. ( 38 ) , and taking the limit as cp -* 0 . The results are: 

I (0) 

r (0) = — ^ ^ la ) 

P C 2 (0) M p (0) 

r p '(o) = { - f(o) H p (o) ^(0) r p (o) + ^ (o) h '( o ) r p (o) 

- c 2 (o) M p (o) r p (o) + ^( 0 ) M p (o) r p (o) 

- l'(0) }/{c 2 (0) H p (0) + c 2 (0) M p (0) j . (to) 


In Eqs . (37) and ( 38 ), the quantities H , M , and I p are functions 
of the steady- state flow variables in the nozzle and these must be com- 
puted before performing the numerical integration to obtain and r • 
For a specified nozzle profile, the steady-state quantities are computed 
by solving the quasi -one -dimensional isentropic steady-state equations 
for nozzle flow. Figure 2 shows the nozzle profile used in our compu- 
tations. All of the length variables have been non-dimensionalized with 
respect to the radius of the combustion chamber, to which the nozzle is 
attached, and hence r = 1. At the throat r.. is fixed by the Mach num- 
ber at the nozzle entrance plane. Hie nozzle profile is smooth and is 




Figure 2. Nozzle Profile Used in Calculating Admittances.- 


completely specified "by r , r and 9 n > which are respectively the radius 
of curvature at the chamber, radius of curvature at the throat and slope 
of the central conical section. The steady-state equations are inte- 
grated using equal steps in steady-state potential <p by beginning at the 
throat and continuing to the nozzle entrance where the radius of the wall 
equals 1 . 

Computations of the admittance coefficients have been performed using 
a three-term series expansion consisting of the first tangential, second 
tangential and first radial inodes . An Adam-Bashforth predictor-corrector 
scheme was used to perform the numerical integration, while the starting 
values needed to apply this method were obtained using a fourth order 
Runge-Kutta integration scheme. The integration computer program has 
been written so that the integration can be performed up to the nozzle 
entrance and also inside the combustion chamber for any desired distance. 
Thus, the admittance relation is obtained at the nozzle entrance section 
or at any station inside the chamber. Computations have been performed for 
several nozzles, at different frequencies and pressure amplitudes of the 
first tangential, mode. 

Figures 3 and 4 show the frequency dependence of the linear ad- 
mittance coefficients for the IT, 2T, and IK modes for a typical nozzle 

(9 s 20°, r = 1.0, r . - 0.9234; M = 0.2) . Here, U) is the frequency 
% p cc cx 

of the IT mode, while the frequency of the 2T and 1R modes is 2u> due to 
nonlinear coupling. Hence the real parts of the linear admittance coef- 
ficients for the 2T and IK modes attain their peal: values at a higher 
frequency than that for the IT mode. The linear admittance coefficients 
for the IT mode are in complete agreement with those calculated previously 

14 / \ 

by Bell and Zinn as ejected from Eq. (22) . 

The frequency dependence of the nonlinear admittance coefficient 
for the 2T mode is plotted in Fig. 5 with pressure amplitude of the IT 
mode as a parameter. While the behavior of the linear admittance co- 
efficient depends only upon the frequency of oscillations , the behavior 
of the nonlinear admittance coefficient is seen to depend on the ampli- 
tude of the IT mode. This result is expected, since in Eq. ( 38 ) , I is 

Jl 

a function of the amplitude of the IT mode. As expected the absolute 
values of both F^ and increase with increasing pressure amplitude of 
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the IT mode, 'which acts as a driving force. It is observed that the 

absolute -values of r . and I\ vary similarly -with frequency as the ab- 

r i 

solute values of and . The frequency dependence of the nonlinear 
admittance coefficient for the 1R mode is plotted in Fig, 6 -with pressure 
amplitude of the IT mode as a parameter. 

Figures 7 and 8 show the effect of pressure amplitude upon the 
magnitude of the 2 ’atio of nonlinear admittance coefficient to the linear 
admittance coefficient for the 2T and 1R modes respectively. These re- 
sults clearly indicate that the nonlinear contribution to the nozzle 
admittance is significant and should be included in nonlinear combustion 
stability analyses . 






Imaginary Part of Linear Ac 




Frequency, <ju 


Figure 5* Nonlinear Admittances for the 2T Mode 















23 


REFERENCES 


1. Tsien, H. S. , "The Transfer Functions of Rocket Nozzles, 

Amer ican Rocket Society Journal , Vol. 22, 1952, pp. 139-143* 

2. Crocco, L. and Cheng, 3, 1., Theory of Combustion Instability in 
Liquid Propellant Rocket Motors , Appendix B. AGARD Monograph No. 8, 
Butterv/orths, London, 1953. 

3. Crocco, L. and Sirignano, W. A., "Behavior of Supercritical Nozzles 
Under Three Dimensional Oscillatory Conditions, ,r Princeton University, 
Department of Aerospace and Mechanical Sciences, Report No. 790, 

April 1967 * 

4. Zinn, B. T. , "A Theoretical Study of Nonlinear Transverse Combustion 
Instability in Liquid Propellant Rocket Motors," Princeton University 
Department of Aerospace and Mechanical Sciences, Report No, 732, 

May 1966. 

5. Zinn, B. T. and Crocco, L. , "Periodic Finite -Amplitude Oscillations 
in Slowly Converging Nozzles," Astronautica Acta , Vol. 13, 1968, 
pp. 481-488. 

6. Zinn, B. T. and Crocco, L. , "The Nozzle Boundary Condition in the 
Nonlinear Rocket Instability Problem," Astronautica Acta , Vol. 13 , 1968, 
pp . 489-496 . 

7. Lores, M. E. and Zinn, B. T., "The Prediction of Nonlinear longitudinal 
Combustion Instability in Liquid Propellant Rockets," NASA CR- 120904, 
April 1972. 

8. Lores, M. E. and Zinn, B. T . , "Nonlinear Longitudinal Combustion 
Instability in Rocket Motors," presented at the AIAA 11th Aero- 
space Sciences Meeting, January 1973* 

9 . Zinn, B. T. and Powell, E. A., "Nonlinear Combustion Instability in 
Liquid Propellant Rocket Engines," Proceedings of the 13th Symposium 
^international) on Combustion , The Combustion institute, pp. 491-503. 

10. Powell, E. A. and Zinn, B. T., "The Prediction of the Nonlinear 
Behavior of Unstable Liquid Rockets," MSA CR-72902, July 1971* 

11. Powell, E. A. and Zinn, B. T., "The Prediction of Nonlinear Three- 
Dimensional Combustion Instability in Liquid Rockets with Conventional 
Nozzles," NASA CR-121279, October 1973* 



29 


‘3i2r- • F inlays oh, :B.vA. and.' S ©riven* L. E., "The Method 1 of Weighted -v •* 
Residuals — A Review," Applied Mechanics Reviews , Vol. 19, Ho. 9, 
September 1966' } pp. 735 “7^9- • ” r ~" T ~~~ 

13. Ames, W. F . , nonlinear Partial Differential Equations in Engineering , 
Academic Pr e s s~Pllew Y ork , 196 5 , PP • ~^*3~262 . 

Ik. Bell, W. A. and Zinn, B. T., "The Prediction of Three-Dimensional 
Liquid-Propellant Rochet Nozzle Admittances," NASA CR- I P II P 9, 

February 1973- 


